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The structure and decay dynamics for 0'^ and 2'*' continuum excitations of are investigated in 
a three-body ^'^O-l-n-l-n model. The validity of a simple approximation for the cross section profile 
for long-lived 2n emission is demonstrated. A sequence of three 0*^ monopole (“breathing mode” 
type) excited states is predicted. These states could probably be interpreted as analogues of Efimov 
states pushed in the continuum due to insufficient binding. The calculated energies of the 2^ states 
are related to the excitation spectrum of We discuss the correlation between the predicted 
spectrum and experimental observations. 
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I. INTRODUCTION 

The interest in extremely heavy oxygen isotopes is 
very high today Ref. EMI- One of the strongest mo¬ 
tivations emanates from the structure theory [Q-Q. The 
evolution of binding energies along the oxygen isotopic 
chain appears to be highly sensitive to the details of in¬ 
teractions and procedures used in the modern structure 
approaches thus providing stringent test for their qual¬ 
ity. This field of research concentrates on the short dis¬ 
tances and short-range correlations. Recently it has been 
demonstrated that another source of an inspiration here 
can be connected to the continuum properties of these 
systems and, correspondingly, to the long-range correla¬ 
tions Ref. 0, [TMl ■ 

The question of prospects to search for neutron ra¬ 
dioactivity has been considered in Ref. 0- In contrast 
to at the proton dripline, characterized by large Coulomb 
barriers, in the vicinity of the neutron dripline the emis¬ 
sion of neutrons from the particle-unstable ground state 
(g.s.) can be hindered mainly by centrifugal barriers (we 
do not consider the possible structural hindrance fac¬ 
tors, connected with many-body effects for the ground 
state neutron emitters). For isotopes with odd num¬ 
ber of neutrons such barriers appear to be sufficient to 
produce long-lived (radioactivity lifetime scale) ground 
states only beyond the s-d shell. Discussion of such 
a remote region of the dripline is not of a practical 
importance nowadays. However, a more complicated 
phenomenon comes into play for the particle-unstable 
dripline systems with even number of neutrons. The pair¬ 
ing interaction can produce specific decay energy con¬ 
ditions which force a simultaneous emission of two (or 
even four) neutrons. The level schemes for and 
illustrating this situation are given in Fig. [T] These, so- 
called true three-body (five-body) decays, are affected 
by the hindrance factor connected with an appearance of 
specific additional barriers in the few-body dynamics [see 
the discussion around Eq. (|31) below]. Thus the search for 
novel types of a radioactivity phenomena, namely, two- 



FIG. 1. (Color online) Experimentally observed energy levels 
of [3 t 3| and Eor theory results the level scheme pre¬ 
dicted for is shown which is based on the one assumed 
for 2®0. 


and four-neutron radioactivities, becomes prospective. 

Among the realistic candidates for 2n and 4n ra¬ 
dioactivity search and were considered in Ref. 
0. Since that time a hint for the very long lifetime 
(T’i /2 = P®) obtained for in Ref. Q- Im¬ 

proved theoretical lifetime estimates in Ref. 0 indicated 
that an extremely low three-body (total) decay energy, 
Et < 1 keV, is required to produce such a long lifetime. 
A realization of such a small decay energy in the nature 
seems unrealistic, however, the experimental decay en¬ 
ergy Et limit for g.s. was steadily decreasing towards 
zero in recent years. At MSU the value Et = 
keV was obtained 0. At GSI Et < 120 keV or Et < 40 
keV was found depending on confidence level 95% or 68%, 
which was chosen for the analysis i- In RIKEN, 26 O and 
2^0 were also produced from 2'^F and 26 f beams with 
much higher statistics than in Refs. Their prelim¬ 

inary data also showed the ground state just above the 
2^0-|-2n threshold, in addition to the newly observed ex- 
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cited state at little over 1 MeV. It should be noted that 
very small decay enerey {Ex ^ 20 keV) was predicted 
theoretically in Ref. [^. Thus, really extreme low decay 
energy of this nucleus cannot be excluded. 

In recent years there were important advances in the 
studies of the lightest 2p emitter ®Be [l3 - [l^ . Consider¬ 
ing also the long history of studies of the 2n halo system 
®He in a three-cluster a+n+n approximation |19l - l^ . it 
can be concluded that the understanding of the three- 
body dynamics in p-shell nuclei is reasonably developed. 
In contrast, for s-d shell systems the situation is much 
less advanced. Only very recently interesting results were 
obtained for the two-proton emitter ^®Ne |^. [^. The 
continuum three-body dynamics of the neuron-rich s-d 
shell nuclei also remains poorly investigated with just few 
examples of such studies [ 1 , llll - [l^ , elaborating mainly 
the ground state properties. 

In previous works on [E[ni we concentrated on 
the decay studies of the g.s. in quite schematic approaches 
aiming first of all at qualitative understanding of the un¬ 
derlying physics for long-living neutron emitters. In this 
work we consider the population and decays of the 
0 "^ and 2^ continuum states using more realistic model 
assumptions, in a broad energy range, and with more de¬ 
tails provided. In this work we report several nontrivial 
results concerning the three-body continuum dynamics of 
the s-d shell nuclei by example of ^® 0 , which provides us 
important general insights in this question. Among the 
obtained results are (i) validity of the simple approxima¬ 
tion to the spectrum shape in the case of radioactive 2n 
decay, (ii) existence of extreme peripheral O'*" monopole 
excitations in the low-lying spectrum of 

II. THEORETICAL MODEL 

The model applied to in this work generally follows 
the approach of Ref. [HI- To provide reasonable pre¬ 
dictions concerning the excitation spectrum of the 
improvements concerning the reaction mechanism treat¬ 
ment were implemented, see e.g. For some direct 

reactions the problem of population and decay of three- 
body states can be formulated in terms of the three-body 
inhomogeneous Schrodingier equation 

— T 3 -t- Vri\-n2 T F^ore-ni T F^ore-n 2 T , 

with the source function $q depending only on one pa¬ 
rameter connected to the reaction mechanism: the trans¬ 
ferred momentum q. The dynamics of the three-body 
^^ 0 -|-n-|-n continuum of is described by the wave 
function (WF) with pure outgoing asymptotic in 

the hyperspherical harmonics (HH) formalism: 

Kj 

XKjip) '^Ki 3 / 2 (^p) ^ exp(-tixp), ( 1 ) 


where H denote the Riccati-Bessel functions of half¬ 
integer index and the “multi-index” 7 denotes the com¬ 
plete set of quantum numbers except for the principal 
quantum number AT: 7 = {L, S, lx, ly}- 

The three-body calculations in the HH method uti¬ 
lize the transition from the three-body Jakobi coordi¬ 
nates {x, y} = {x, fix, y, fly} to the collective coordinates 
{p, Hp} = {p,0p,flx,fly}. The hyperradius p (describing 
collective radial motion) and the hyperangle 9p (respon¬ 
sible for geometry of the system at given p) are defined 
via the cluster coordinates as: 


X = 


f "4iri 

V 


_ , (Ai-|-.42).43 

•y V Ai-l-.42-l-.43 

= , 


+A2r2 


- ra 


+A 2 

9p = aTctaii{x/y). ( 2 ) 


Hypermomentum x = y/2MWr is the dynamic variable 
conjugated to hyperradius. The mass M is an average 
nucleon mass for the considered nucleus. 

The hyperspherical harmonics with definite to¬ 

tal angular momentum J and its projection Mj 


[[FiJH,) 0 Yi^{fly)]L 0 Xs\j^^ 


form a full set of orthogonal functions on the five¬ 
dimensional “hypersphere” flp. The pure hyperangular 
functions are expressed in terms of Jacobi polino- 
mials. 

The three-body Schrodinger equation ([T|) in the hyper¬ 
spherical basis is reduced to the set of coupled differential 
equations for the functions : 


C{C + l) 


dp^ 


+ 2M{E-VKxK-t{p)} 


X^kI ip) 


— 2M 'y ( ^K'f,K''f'{p)x\i'^-f'ip) ‘d'M ^ci^K-fip), (3) 


which can be interpreted as motion of a single “effective” 
particle in a strongly deformed field. The “three-body 
potentials” (matrix elements of the pairwise potentials) 
VK'^,K''y'ip) and the partial source terms are de¬ 

fined as 

VKxK'Yip) = Jdflp J^^/*(Hp) E F„ (ri, ) > 

i<j 

d>q.iC-i(p) = I dflp J^f^riflp) $q(p, Hp) . 


The details of the hyperspherical method application to 
various three-body systems in different ph ysic al situa¬ 
tions can be found in the papers [H, [H, [26l4^ . 

The important qualitative difference between Eq. ® 
and the conventional two-body situation is that the “ef¬ 
fective angular momentum” C — K -\- 3/2 in the three- 
body problem is not equal to zero even for the lowest 
possible quantum state with K = 0. Thus, there exists a 
three-body centrifugal barrier even for decays via s-wave 
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emission of neutral particles producing strong hindrance 
factors for the widths of such low-energy decays. 

The differential cross section is expressed via the flux j 
induced by the WF on the remote five-dimensional 
surface flp with p = pmax 


da 


■J, j = 


,(+)\ 


dEjE dVtp 
1 

= TT I™ 
M 


Pma 




(4) 


pma 


The approach with inhomogeneous Schrodinger equa¬ 
tion o had previously been applied to two differ¬ 
ent direct reaction mechanisms (knockout and charge- 
exch ange) p opulat ing the three-body continuum of the 

®Be [iMil, i°He Q, and ^^Ne [13,HI]. 

The source function <i>q for the 0"^ continuum was ap¬ 
proximated assuming a sudden removal of a d-wave pro¬ 
ton from 


$(o+) = 

q 




(d>24o |d>27F) 


(5) 


where is the radius-vector of the removed proton. The 
^^F g.s. WF 4['27 f was obtained in a three-body 
cluster model and the technicalities of proton removal 
from the ^®F core of the ^^F nucleus are the same as in 
calculations of Ref. . 

The source function <i>q for the 2“*' continuum can not 
be easily evaluated by a simple proton removal model 
in the framework of three-body approach to structure of 
^^F. Therefore, we use the source generated by addition¬ 
ally acting on the valence neutrons of the ^^F g.s. WF by 
the quadrupole operator: 

=V2 f d^r„e*‘’''"(^'24o| ^ r'fY2mi{fi)\'^27p). 

1^1,2 

(6) 

Expressions of this type typically arise in the direct re¬ 
action studies and seem to be sufficiently suited for ex¬ 
ploratory studies of 2+ excitations in the case. The 
a^roximation is the same as used in the recent paper 

The sudden removal approximation is not intended 
for absolute cross section calculations, and therefore the 
source strength coefficients Vi in Eq. and (ID) are ar¬ 
bitrary values providing the source function the correct 
dimension of [energy/length®/^]. 


III. POTENTIALS 

In this work we employ for the nucleon-nucleon chan¬ 
nel the quasirealistic potential from Ref. including 
central, spin-orbit, tensor, and parity-splitting terms. 

In the ^^0-n channel we used the potential from m 
characterized as “moderate repulsion in s and p waves”. 
This is a Woods-Saxon potential with radius rg = 3.5 fm, 
diffuseness a = 0.75 fm, and depth parameters Vs = 70 


> 

<D 






FIG. 2. (Color online) The 0"^ excitation spectrum of 
calculated by the three-body model for proton knockout from 
Weights, Wi, of the peaks are shown in percent. 


MeV and I^ = 70 MeV for s and p waves respectively. 
The d-waves component was modified by an inclusion 
of the Is interaction to produce realistic excitation spec¬ 
trum of Vd = —33 MeV, Vis = —5 MeV. We also 
varied the Is interaction to check the sensitivity to this 
parameter, see Sections IVll and IVlIl 

The three-body potential V 3 (p) depending on p with 
the Woods-Saxon parameterization (po = 5 fm, a = 0.9 
fm) is used to control the decay energy Et when the fine 
adjustment is needed, see the discussion in [3l|. This 
potential is very small (on the level of just a few keV) for 
the potential chosen in this work and therefore V^ 

does not influence the other results of the calculations on 
the practical level. 


IV. EXCITATION SPECTRUM AND 
STRUCTURE OF 0+ STATES 


A calculated spectrum for 0+ continuum states of 
is shown in Fig. [2] In this particular calculation the g.s. is 
obtained at Et = 10 keV. Two more states are obtained 
at about 1.76 and 2.6 MeV. Some hint for a broad 0+ 
state at about 4.2 MeV can also be seen. Within the sim¬ 
ple reaction model used in this work the ground 0+ state 
is populated with the relative intensity Wi ^ 92% indi¬ 
cated in Fig. ID while the excited 0+ states are populated 
at a level of Wi>i ~ 2 — 3 %. It can be expected that 
in more complicated reaction scenario the source func¬ 
tion <I>q should have smaller affinity to the ground 
state WF and we can expect up to Wi>i ^ 5—15 % 
population rates for the excited 0+ states. The rela¬ 
tive population intensities above are calculated by in¬ 
tegration of the strength function Fig. [5] in the range 
[Esrii) — r(i), Esrii) + r(i)] around the energy E^rii) of 
the i-th resonance which width is F(i). 
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V. R-MATRIX LIKE PHENOMENOLOGY 


The ideas of using R-matrix type expressions for anal¬ 
ysis of the three-body excitations and decays have been 
discussed occasionally in the literature [ij, [H, in 
quite a sketchy way. It is important to understand gen¬ 
eral validity and limits of applicability of such expressions 
for practical application as for theoretical estimates and 
as for phenomenological analysis of experimental data. 

It can be shown that for the source function normalized 
at given q value 

J dpp^dnp2M\^^{p,nX = 


the energy proHle of flux is provided by the expression 

TT FkIEt) 

^ {Et - E^r? - T{E:ir?ld ’ 
2^ _ (2/7r) 

^ ^K+2\^h'cii) ^ ^ K+2 


j{Et) = W{E3r) 


^k{Et) = 2 7wl -2— - . 2 —1 -T ■ 

K+2\^P<^'^) E ^^4-2 (-^Pch) 


(7) 


Here the hypermomentum x — y/2ME^ is defined at 
the three-body resonance energy E^r, and the “Wigner 
limit” 7wL = estimates the upper limit for 

the width. The functions J and Y are cylindrical Bessel 
functions regular and irregular at the origin, respectively. 
This expression is totally analogous to the standard two- 
body R-matrix expression as it is based on the assump¬ 
tion that the penetration in the three-body system is 
defined by the three-body channel in the hyperspherical 
decomposition of the WF with the lowest possible hyper¬ 
momentum R = Kmin- For systems with zero spin of the 
core the latter is trivially related to the total spin of the 
state: iFmin = J. The “normalization” W{E^r) is con¬ 
nected to the affinity of the source to the inner structure 
of the resonance; this value is expected to be smaller than 
unity, but of the order of unity for the realistic situation. 
It was found to be IF(g.s.) = 0.92 in our calculations, see 
also Fig. [5] 

The example of application of the R-matrix expression 
is provided in Fig.|31 The possibility of fitting the excita¬ 
tion profile by the Breit-Wigner shape in the vicinity of 
the resonance does not cause any doubts. However, it can 
be seen that the near-perfect description is provided by 


TABLE I. The three-body structure of the g.s. of initial 
nucleus and 0^ excitations in terms of probability 
in percent of the corresponding [Z|]o configuration. The en¬ 
ergies of the states with respect to the three-body breakup 
threshold are provided in the last row. 


Ij 

'"F, g.s. 

'®0, g.s. 

0+ 

2^0, 0+ 

Sl/2 

0.55 

0.67 

3.7 

3.8 

^3/2 

84 

79 

80 

86 

^5/2 

13 

19 

6.0 

6.1 

Et (MeV) 

-3.2 

0.01 

1.7 

2.6 



FIG. 3. Cross section profiles for the very narrow g.s. 
calculated by the three-body model and obtained by the R- 
matrix-type approximation are shown on two different scales 
in panels (a) and (b). 


the R-matrix expression up to something like 10® widths 
away from the resonance. For this Figure the “spectro¬ 
scopic factor ” is fitted in such a way that the cal¬ 
culated width Tk for the three-body resonance energy 
E^r = 9.7 keV is exactly reproduced by the R-matrix 
expression [in this specific case TK[E^r) = 4.04 x 10“® 
keV]. 

It is also necessary to fix one more parameter: the 
“channel radius” pch- In the three-body case this pa¬ 
rameter does not have such a well defined meaning as 
in two-body R-matrix phenomenology and some investi¬ 
gation is required here. Definition of the spectroscopic 
factor for the single-channel approximation is 

0\^{pc\,) = = ^ j^dx\xK'^{xpc\,)?. 


The dependence on channel radius, obtained using 
the calculated three-body WF, is illustrated in Fig. 
m It is very stable in a broad range of channel radii, 
varying from 0.007 to 0.004 for 4 < pch < 25 fm. The 
6^ variation for the mentioned range of channel radius is 
in very good agreement with the relative weight of the 
[s^] configuration calculated within the whole internal 
region VF(s^) = 0.0067, see also Table ID The calculated 
three-body width is reproduced for pch = 13 fm and for 
less than ±50% variation of width we need to keep the 
channel radius in the range 9 < Pch < 16 fm. 

We can see that the extension of a simple single chan¬ 
nel R-matrix phenomenology, Eq. 0, to the three-body 
decays provides easily tractable and very reliable results 
for long-lived two-neutron emitters. We can also con¬ 
clude that the use of typical R-matrix parameters, cho¬ 
sen according to the prescription discussed above, can be 
expected to provide widths values with an uncertainty 
around 50%, which is a quite accurate result for esti¬ 
mates concerning long-living states. 
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FIG. 4. (Color online) Spectroscopic factor 6^ Eq. ((7|) as a 
function of channel radius pch- The dotted curve shows the 
dependence of the “internal normalization” I on pch- 


VI. THE MONOPOLE 0+ EXCITATIONS 



FIG. 5. (Color online) WF densities for the ^®0 ground state 
and two excited O’*" states. All WFs are normalized to unity 
maximum value. Arrows indicate the humps connected with 
monopole excitations. 


The nature of the excited 0+ states is very interesting 
and deserves a special discussion. Table H] provides the 
basic structure information about the O'*" states in¬ 
dicating their high similarity. It is clear that a simple 
explanation for such a situation is that the predicted ex¬ 
cited 0+ states are all monopole (often called “breathing 
mode”) excitations. 

To check this assumption we have studied the radial 
evolution and the correlation densities of the O'*" 
WFs at corresponding energies. One can see in Fig. [5] 
that the g.s. WF density decreases more or less exponen¬ 
tially inside the barrier. At larger distances the behav¬ 
ior tends to be constant which corresponds to approach¬ 
ing the asymptotic behavior ^ exp{+i>cp) of the 

three-body WF. However, for the Oj and state WFs 
there exists one and two extra humps respectively (indi¬ 
cated by arrows in Fig. ED in the above-the-barrier slope 
of the density before the asymptotic (constant) behavior 
is achieved. This is exactly expected for the monopole 
states where the major WF component should have one 
or more nodes in the radial WF. The extreme radial ex¬ 
tent of the humps is notable. E.g. the peak for Oj with 
p ^ 30 fm corresponds to typical single particle distances 
in the core-neutron channel of about 20 fm. 

To understand the reasons for the monopole state for¬ 
mation we performed simple width estimates for major 
configurations of the WF. Fig. E] provides the esti¬ 
mates in a “direct decay” R-matrix model from Ref. [3l| : 


Jl32 


{Et) = 


Et {Vzf 

27r 


de 


Tj^^eEr) 


{eET-E,,Y+T,,{eETYlA 


T,,{{l-e)ET) 


((1 _ e)ET - E,,Y + r,,((l - e)ETYlA ' 


( 8 ) 


This model is reasonably well approximating the true 
three-body decay mechanism, and also provides a smooth 
transition to the sequential decay regime. The direct 


decay model is constructed in the spirit of independent 
particle approach with two nucleons being emitted from 
states with definite single-particle angular momenta ji 
and j 2 sharing the total decay energy Et, but with in¬ 
teraction between nucleons neglected. The F^-^ is the 
standard R-matrix expression for the width as a function 
of the energy for the involved resonances in the ^"'^0-l-n 
subsystems. It is assumed that two-body resonant states 
with energies Ej. are present in both two-body subsys¬ 
tems and that the values F^. [Ej. ) correctly describe their 
empirical widths. The matrix element (V3) can be well 
approximated as 


(V3)" = D:,[{ET-E,,-E,,Y + {T,,{E,,)+T,,{E,,)flA\ , 


where the parameter ZI3 k. 1.0 — 1.5 is a constant. 

It can be seen in Fig. El that for low-energy decays of 
the QT states the predominant contribution to the width 
is connected with the configuration of Here 

the width associated with the dominant [d3/2]o conhg- 
uration of the WF (see Table HD is strongly suppressed. 
However, for Et > 0.77 MeV the decay mechanism for 
the [d3/2]o configuration changes to the sequential emis¬ 
sion via the ^3/2 ground state of So, at higher ener¬ 
gies the partial width for this configuration rapidly grows 
and becomes approximately equal to that of the [s^^2]o 
conhguration at Et = 1.6 MeV. This is approximately 
the energy at which the second 0+ state appears. So, 
despite the simplicity of these estimates they provide a 
hint that the appearance of the excited QT states could 
be connected with the ability of the [^3^210 configuration 
to propagate effectively to large distances above the bar¬ 
rier. We can estimate that the monopole states found, are 
built on interference patterns between [s3^2]o and [d3^2]o 
components at large distances. 

To get a deeper insight into the process we may inves- 
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Ej (MeV) 

FIG. 6. (Color online) Width estimates for the O’*" and 2'*’ 
states of in the R-matrix-type approach of Eq. ([HI) for 
different [ 3132 ] j configurations. Vertical dashed lines indicate 
the predicted positions of the 0^ states. 


tigate the three-body WF correlation densities Fig. [7J 

w{p, 0 p) = J dQy\^‘-+\p,9p,n^,ny)\'^. 

In the “Y” Jacobi system we should choose the mass 
number of the cluster as Ai or A 2 in Eq. ([2]) and for 
such a relatively heavy core cluster we can approximate 
the distances between core and valence nucleons by: 

^core-ni = \JP sin(6»p), 

?-core-„. ~ P COs(0p). 

Thus, the provided correlation densities illustrate the 
evolution of relative distances in the single-particle chan¬ 
nels. 

It can be seen in Fig.[7]that the correlation densities for 
the calculated Oj and states have a complicated cor¬ 
relation pattern in the hyperangle Op at large distances 
10 < p < 30 fm. Such a triple-peak pattern should be 
connected with the important contribution of the [d^] 
configuration in the above-the-barrier region. The 
configuration is as expected dominant in the nuclear in¬ 
terior p < 10 fm, but for the ground state it is suppressed 
under the barrier, see Fig. [3 (a). In contrast, for the 0^ 
and O3 states, illustrated in Figs. [3(b) and (c), the [d^] 
component extends also to the peripheral region forming 
the triple-peak patterns. 

Figure |5| demonstrates the evolution of the 0"^ states 
with Vis- This evolution provides some confirmation that 
the “remnants” of the broad monopole O4 state is present 
in the spectrum Fig. |3 at about 4 MeV. With Vis ap¬ 
proaching zero, the O'*" spectrum shifts downwards in en¬ 
ergy providing a much narrower, and thus well defined, 
O4 state. For larger negative Vs values the O4 state 
increases in energy and is totally dissolved in the contin¬ 
uum. 
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FIG. 7. (Color online) The WF correlation density 
W{p,0p) in the “Y” Jakobi system for the ground state (a) 
and two excited O'*" states (b) and (c). Notice the change of 
p scale. 


VII. THE 2+ STATES 

The excitation spectrum of is expected to have 
relatively low level density due to the simplicity of the 
spectrum of ^^0, where only one low-lying state is known 
so far (d3/2 at 0.7 MeV). In the calculations we obtain 
just two 2+ states where the lower one has a structure 
characterized by the configuration mixing [d3^2]2"['^5/2]2j 
while the higher one is based on the {si/ 2 d^p^ 2 -{si/ 2 db/ 2]2 
configurations. Within the three-body model the sepa¬ 
ration between these configurations is defined by the Is 
splitting between the ^^0 g.s. d ^/2 and the yet unknown 
d ^/2 state. Figure |5| demonstrates the dependence of the 
positions of the 2+ states on the intensity of the Is inter¬ 
action under the condition that the d ^/2 g-S. position 
remains fixed. The 2+ states naturally become degener¬ 
ate for Vis = 0. 

Most of the calculations of this work are performed for 
Vis = —5 MeV, which provides the position of the 0+ g.s. 
to be exactly on the 2n threshold. Under this condition 
the predicted positions of the 2+ state are 1.6 MeV and 
4.5 MeV. The calculated three-body width value for the 
2^ state is F = 115 keV. However, the width of the states 
with expected dominant sequential decay mechanism are 
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FIG. 8. (Color online) The systematics of the 0^ and 2'*’ 
states of as a function of the Is interaction intensity Vis 
in the d-wave ^“^O-n channel nnder the condition that the 
^ 3/2 g-S. position remains fixed. The result for Vis = —5 MeV 
(indicated by vertical dots) is provided in Fig. [T]as our main 
result. The predictions of Ref. are indicated by arrows in 
the left part of the plot. 

not reliably predicted in the HH method calculations. 
For that reason we also performed estimates using the 
direct decay model expression Eq. The result for the 
[si/ 2 d 3 / 2]2 configuration which has the most favorable 
penetration conditions is provided in Fig. [51 For the 2^ 
state the value Fi/ 2 . 3/2 ~ 2.7 MeV. However, it should be 
corrected for the weight of the [ 51 / 2 ^ 3 / 2)2 configuration 
in the interior of the three-body WF 

r = ■ 

The value VFi/ 2 , 3/2 = 3.6%, which is quite small, is ob¬ 
tained in the three-body calculations. Then the estimate 
F = 96 keV is obtained. This value agrees well with the 
results of the three-body width calculations, and we es¬ 
timate that it is realistic to expect the width of the first 
2+ state to be in the range 100 — 120 keV. 


VIII. DISCUSSION OF THEORETICAL 
RESULTS ON ^^O SPECTRUM 

In our calculations we have chosen the calculation 
scheme, where the effect of the occupied ^ 5/2 orbital in 
is imitated by inverse Is forces moving the d ^/2 state 
to an energy higher than that of the ^ 3/2 state. The par¬ 
ticular value of this interaction was fixed by reproducing 
the near-threshold position of the g.s. Such a com¬ 
putation scheme is not free of problems, but it seems, 
however, to be quite successful for the nucleus. This 
nucleus is just two neutrons away from the ^®0 neutron 
shell closure and the main dynamical degrees of freedom 
are connected just with the motion of the two valence 
nucleons. This is a clear motivation for the use of the 
three-body core-|-n-|-n model. 


The oxygen isotope chain was studied in details by 
Volya and Zelevinskiy using the continuum shell 
model. They predicted the ground state position to 
be practically exactly on the threshold with Et =21 
keV. Such a bold prediction was nicely confirmed by 
the recent experimental studies S i 0. The predic¬ 
tions of @ about 2+ states are in reasonable agreement 
with our calculations: our 2 ^ state position is around 
100 keV lower in excitation energy than in Ref. Q and, 
correspondingly, the 2^ is about 1.7 MeV lower. The 
^i>i low-lying excited states are absent in the contin¬ 
uum shell model calculations of Ref. @. This, probably, 
could be connected to the “short range character” of the 
shell model calculations in general: in our calculations 
the excited O’*" states require the dynamical range of tens 
of Fermi to be properly accounted for. Another possible 
reason could be that the treatment of ^'^0+n+n con¬ 
tinuum in Ref. @ is not fully dynamical, as it relies on 
the simplified three-body Green’s functions not including 
N-N final-state interaction. For this reason it is a aston¬ 
ishing that the width obtained in the current calculations 
practically coincides with results the of @ . According to 
our experience, the simplified three-body Green’s func¬ 
tion should provide a width for the [s^] configuration de¬ 
cay which is ^ 10 — 30 times smaller than that obtained 
if appropriately taking the final state nucleon-nucleon in¬ 
teraction into account. 

The g.s. was studied in [T^. The higher 0+ exci¬ 
tations are not discussed in that work but there are some 
indications of the excited 0+ at about 3.4 MeV, see. Fig. 
1 of Ref. [l2|. The authors do not elaborate this result 
and its reliability is not known; at least there is some in¬ 
dication of the possibility of an existence of the low-lying 
excited in an alternative approach. The first 2+ state 
was obtained at about Et = 1.35 MeV in the recent cal¬ 
culations Ref. M- This value is not drastically different 
from our results and the results of 

The ab-initio shell model theoretical calculations of 
Ref. [ 3 ^ provided the excitation energy of the 2+ state 
in in the range E* = 1.2 — 1.7 MeV depending on 
the details of calculations. The obtained results are rea¬ 
sonably consistent with other theoretical predictions. 


IX. DISCUSSION OF RELEVANCE TO 
EXPERIMENTAL DATA 

How could the predicted excitations of in our work 
be related to observations? The considerable evolution of 
the spectrum shape with the Vis parameter is shown 
in Figure jS] However, it should be noted that the pre¬ 
dicted picture of excitation energies is relatively stable 
in this plot for reasonable variation of Vs- In particular, 
the excitation energies for states and 2 ^ state are 
practically constant. 

There is some evidence for intensity at about Et = 
2 MeV in Refs. S 0. A peak just above 1 MeV was 
observed in Ref. [ 3 . For this energy range we predict 







relatively narrow (F ~ 0.11 MeV) 2^ state at Et ^ 1-6 
MeV. The right “wing” of this peak could be situated on 
a “background” formed by 0 ^ and, maybe, states. 

Evidence for the excited state at about Et = 4.2 
MeV was obtained in Ref. Q, although with marginal 
statistics. Within this energy range we predict quite 
broad and overlapping 2 ^ and O 4 states. 


X. CONCLUSIONS 

We studied the 0+ and 2+ continuum properties of the 
system in three-cluster ^^ 0 -l-n-l-n theoretical model. 
The main results obtained in this work are: 

(i) The ground state decay spectrum of the long-living 
2 n emitters can be very well approximated in a broad en¬ 
ergy range by a simple analytical expression generalizing 
R-matrix phenomenology for use in the hyperspherical 
space of three-body systems. 

(ii) A number of monopole (breathing mode) excited 0^ 


states are predicted in at about 1.76, 2.6, and 4.2 
MeV. These states extend to extreme distances in radial 
space with typical single-particle orbital sizes around 20 
fm. The predicted densely spaced sequence of 0-1- states 
in the continuum states with such radial proper¬ 
ties resembles the Efimov phenomenon. These states can 
probably be considered as the lowest Efimov states forced 
into continuum by insufficient binding. 

(iii) The predicted 2+ states in at 1.6 and 4.5 
MeV are in a reasonable agreement with continuum shell 
model calculations. 

(iv) We suggest to identify the experimentally observed 
intensity in the range 1 — 2 MeV in the spectrum of as 
a “pileup” of Oj and 2f state contributions. The possible 
experimental peak at ^ 4 MeV in can be associated 
with overlapping broad 2 j and continuum states. 
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